NASA Technical Memorandum 109150 



//v -;8 

&Z2 2 - 


Three-Axis Stabilization of Spacecraft 
Using Parameter-Independent Nonlinear 
Quaternion Feedback 


Suresh M. Joshi 

Langley Research Center , Hampton, Virginia 
Atul G. Kelkar 

NRC Fellow, Langley Research Center, Hampton, Virginia 


( NAS A- TM- 109 150) THREE- AX I S N95-1 3729 

STABILIZATION OF SPACECRAFT USING 
PARAMETFR-INDEPENOENT NONLINEAR 

QUATERNION FEEDBACK (NASA. Langley Unclas 

Research Center) 10 p 


September 1994 


G3/18 0028222 


National Aeronautics and 
Space Administration 
Langley Research Center 
Hampton, Virginia 23681-0001 



Summary 


This paper considers the problem of three-axis attitude stabilization of rigid spacecraft. 

A nonlinear control law which uses the feedback of the unit quaternion and the measured 
angular velocities is proposed and is shown to provide global asymptotic stability. The 
control law does not require the knowledge of the system parameters, and is therefore robust 
to modeling errors. The significance of the control law is that it can be used for large-angle 
maneuvers with guaranteed stability. 

Introduction 

Attitude control of a free-flying spacecraft has long been known as an important 
problem, and has been the subject of many papers since the fifties and sixties [1,2,3]. It 
is also a unique problem in dynamics because of the fact that the finite rotation of a rigid 
body does not obey the laws of vector addition (in particular, commutativity) and, as 
a result, the angular velocity of the body can not be integrated to give the attitude of 
the body. The most widely used method of defining the rotation of a body between two 
different orientations is an Euler angle description. A 3 x 3 direction cosine matrix (of 
Euler rotations) is used to describe the orientation of the body (achieved by three successive 
rotations) with respect to some fixed frame of reference. However, there is an inherent 
geometric singularity in the Euler representation. This problem can be avoided by using a 4- 
parameter description of the orientation [2,3,4], known as ‘quaternions’, which can be used 
to describe all possible orientations. The quaternion approach uses Euler’s theorem which 
states that any rotation of a rigid body can be described by a single rotation about a fixed 
axis. The advantage of using quaternions is that successive rotations result in successive 
quaternion multiplications which are commutative. Some early results on the use of 
quaternion feedback for attitude error representation and automatic control of the attitude 
can be found in [1]. Quaternions were used for the simulations of the rotational motion of 
rigid bodies as early as the 1950’s [5], The use of quaternion feedback for controlling robotic 
manipulators can be found in [6,7], and for spacecraft control can be found in [8-11], 


Various linear and nonlinear quaternion-based control laws have been recently proposed 
[10, 11] for the attitude control of a single-body rigid spacecraft. However, the control laws 
proposed in [10] require the knowledge of the system’s moments of inertia and also constrain 
the choice of the gain matrices. In [11], both model-dependent and model-independent 
control laws were presented; however, the control laws used scalar gams. 

In this note, a model-independent, nonlinear control law is presented which uses 
quaternion feedback and symmetric and positive definite gain matrices. Global asymptotic 
stability of the proposed control law is shown by using Lyapunov analysis. The Lyapunov 
function used here for proving asymptotic stability does not need a cross term similar to the 
one used in [11], and the proof is made much simpler. 

Quaternion Feedback Control 

The rotational equations of motion of a rigid spacecraft are given by. 

Ju + ux (Jw) = u (1) 

where J is the 3 x 3 inertia matrix; - is the 3 x 1 angular velocity vector; and u is the 3 x 1 
vector of actuator torques. The objective of the control system is to bring the spacecraft to 
the desired attitude (orientation) starting from any initial condition. 

The orientation of a free-floating body can be minimally represented by a 3-dimensional 
orientation vector. However, as stated previously, this representation is not unique. One 
minimal representation that is commonly used to represent the attitude is Euler angles. 

The 3x1 Euler angle vector rj is given by : £(»?)»? = w, where £(»?) is a 3x3 transformation 
matrix. E( V ) becomes singular for certain values of „; however, it should be noted that 
the limitations imposed on the allowable orientations due to this singularity are purely 
mathematical in nature and do not represent physical restrictions. The problem of singularity 
in a 3-parameter representation of attitude has been studied in detail in the literature 
[2,3,8,10,11]. An effective way of overcoming the singularity problem is to use the quaternion 
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formulation (see [1]- [3]). The unit quaternion a is defined as follows. 


a 4 = cos(^) (2) 

a = (d 1 ,d 2 ,Q3) T is the unit vector along the eigen-axis of rotation and <j> is the magnitude of 
rotation. The quaternion is also subjected to the norm constraint: 

a T a + a 2 4 = 1 ( 3 ) 

It can be also shown [12] that the quaternion obeys the following kinematic differential 
equations. 


= i(w x a + a 4 w) 

(4) 

1 T— 

° 4 = ~2 U a 

(5) 


The quaternion can be computed [12] using Euler angle measurements (Eq. 4). The 
equilibrium solutions of the open-loop system, given by equations (1), (4), and (5), can be 
obtained by setting all the derivatives to zero. That is, 


a = {a T ,c * 4 } 3 


ai 
<* 2 
<*3 


us x (Jus) = 0 

(6) 

wxa-f a 4 w = 0 

(7) 

us T a = 0 

(8) 


Taking the dot product of u> with both sides of Eq. (7), 


a 4 (u; • u>) = 0 

That is, either ct 4 = 0, or u — 0, or both are 0. If a 4 ^ 0, then u> = 0. If c* 4 = 0, then from 
(7) and (8), w x a = 0, and a; • a = 0, i.e., lj = 0, or a = 0, or both are zero. However, from 
( 3 ), a 4 = 0 =i> a ^ 0; therefore, w = 0 when the system is in equilibrium. The system has 
multiple equilibrium solutions: (a SS T , a 4ss ), where, the subscript ‘ss’ denotes the constant 
steady-state value. 
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Consider the control law u, given by: 


u = — i[(a + a 4 /)G p + 7(1 - a 4 )f]a — G r w (9) 

where G p and G r are symmetric positive definite (3 x 3) matrices; 7 is a positive scalar; and a 
represents the 3x3 cross product matrix of the vector a. Equation (9) represents a nonlinear 
control law. The following result gives the closed-loop equilibrium solutions. 

Lemma 1. Suppose G p is symmetric and positive definite and 0 < Am(G p ) < 27, where A a/(0 
denotes the largest eigenvalue. Then the closed-loop system given by (1), (4), (5) and (9) 
has exactly two equlibrium solutions: [a = w = 0, a 4 = 1] and [a = w = 0, a 4 = -1]. 

Proof.- The closed-loop system is in equilibrium when the derivatives in equations (1), (4), 
and (5) are zero. Proceeding as in the open-loop case, the closed-loop equilibrium solution is 
given by: w = 0,c? = a ss ,a 4 = a 4ss . From equation (1), ^ = 0 => u = 0. From equation (9), we 
have 

[( 5 +a 4 I)G p + y(l-a 4 )I]a = 0 (10) 

Pre-multiplying the above equation by a T and noting that the first term vanishes, we have 
the following: 

u = 0 => qT Ma = 0 


where 

M = a 4 G p + 7(1 — a 4 )/ (11) 

The eigenvalues of M are given by: A i(M) = a 4 \i(G p ) + 7(1 - ar 4 ) = « 4 (A,(G P ) - 7) + 7- M is 
singular when A, (AT) = 0, i.e., when a 4 = There are three different subcases that 

need to be examined: (a) 0 < A,(G P ) < 7, (b) A,(G P ) = 7, and (c) 7 < A,(G P ) < 27. In subcases 
(a) and (c), A ;(M) = 0 only if |a 4 | > 1, which is not feasible, since -1 < a 4 < 1. That means, 
for subcases (a) and (c), A ,(M) ^ 0 for any feasible values of a 4 - In subcase (b), i.e., when 
Aj(G p ) = 7, A,(Af) = 7 > 0. Therefore, M is nonsingular, and a = 0. Then, from equation ( 3 ), 
we have: a 4 = 1 or -1. (Note that, if AdG p ) > 27 then there are feasible values of a 4 for which 
Ai (AT) = 0). ■ 
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(The symbol ■ denotes the end of the proof.) 

From Lemma 1, there appear to be two closed-loop equilibrium points corresponding to 
a 4 = 1 and a 4 = -1 (all other state variables being zero). However, from equation (2), 
a 4 = 1 => <j> = 0, and a 4 = -1 => <j> = 2tt (or more generally 2mr), i.e., there is only one 
equilibrium point in the physical space. We shall define the desired equilibrium state as: 
a = w = 0 , 0-4 = 1. In order to make the origin of the state space the desired state, define 
0 = (a 4 - 1). Equations (4) and (5) can then be rewritten as: 

a = |(uxtt + (/J+l)u) (12) 

0=-^u> T a (13) 

The system represented by equations (1), (12) and (13) can be expressed in the state-space 
form as follows: 

x = f(x,u) (14) 

where x = (a T ,0, u> T ) T . Note that the dimension of x is 7, which is one more than the 
dimension of the system. However, one constraint (Eq. 3) is now present. It can be easily 
verified from (4) and (5) that the constraint (3) is satisfied for all t > 0 if it is satisfied at 
t = 0. 

If the objective of the control law is to transfer the state of the system from one orientation 
(equilibrium) position to another orientation (i.e., a rest-to-rest maneuver), then without 
loss of generality, the target orientation can be defined to be zero. The initial orientation, 
given by (a(O),0(O)) can always be defined in such a way that — 1 < 0 < 0 (i.e., 0 < a 4 (0) < l), 
corresponding to \<f>\ < n. 

Consider the control law given by: 

u = -^[@+(0+l)I)G p - 7 0I}a-G r u (15) 

The control law given by (15) was stated in [11]; however, conditions for the existence of the 
closed-loop equilibrium solutions were not investigated, and stability proof was not given. 
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The following theorem establishes the global asymptotic stability of the physical 
equilibrium state (the origin of the state-space) of the system. 

Theorem 1. Suppose G v and G r are symmetric and positive definite, and 0 < A A /(G p ) < 

2j. Then, the closed-loop system given by equations (1), (12), (13), and (15) is globally 
asymptotically stable (g.a.s.). 

Proof. Consider the candidate Lyapunov function 

V = uT Jui + a T G p a + y/3 2 (16) 

V is clearly positive definite and radially unbounded with respect to the state vector 
x = {a T ,p,w T } T . Taking the time derivative of V, we have: 

V = 2 w T [-w x ( Ju> ) + u] + 2 a T G p (u x a + {0 + l)w) - jpui T a (17) 

Noting that the first term in V is zero, and substituting for u from (15), after simplification, 
we get: V = -w T G T w, i.e., V is negative semidefinite. V = 0 only when w = 0. Following the 
same procedure as Lemma 1, it can be shown that V - 0 only at the two equilibrium points, 
a = w = 0,/? = 0 (corresponding to o 4 = l) and q = w = 0,/? = — 2 (corresponding to a 4 = -l). 

Consistent with the previous discussion, these values correspond to two equilibrium 
points representing the same physical equilibrium state. It can be easily verified, from 
equation (16), that any small perturbation e in (3 from the equilibrium point corresponding 
to /? = -2 will cause a decrease in the value of V (e has to be > 0 because -2 < 0 < 0). Thus, 
in the mathematical sense, /? = — 2 corresponds to an isolated equilibrium point such that 
V - 0 at that point, and V < 0 in a neighborhood of that point, i.e., /? = -2 is a Tepeller’ and 
not an ‘attractor’. It has been already shown that V is negative everywhere in the feasible 
state space except at the two equilibrium points. That is, if the system’s initial condition 
lies anywhere in the state space except at the equilibrium point corresponding to /? = -2, 
then the system will asymptotically approach the origin (x = 0); and if the system is at 
the equilibrium point corresponding to 0 = — 2 at t = 0 then it will stay there for all t > 0. 
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However, this is the same equilibrium point in the physical space; hence, it can be concluded 
by LaSalle’s theorem that the system is globally asymptotically stable. ■ 

Concluding Remarks 

The problem of three-axis attitude stabilization of a spacecraft was considered. A 
nonlinear quaternion- based feedback control law was given, and was shown to provide 
global asymptotic stability. The control law does not depend on the knowledge of the 
system parameters (i.e., moments of inertia), and is therefore robust to modelling errors and 
parametric uncertainties. 
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